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We calculate the nuclear modification factor of neutral pions and the photon yield at high pr 
in central Au-Au collisions at RHIC (s/s =200 GeV) and Pb-Pb collisions at the LHC (yfs =5500 
GeV). A leading-order accurate treatment of jet energy loss in the medium has been convolved with 
a physical description of the initial spatial distribution of jets and a (1+1) dimensional expansion. 
We reproduce the nuclear modification factor of pion Raa at RHIC, assuming an initial temperature 
T, =370 MeV and a formation time r, =0.26 fm/c, corresponding to dN/dy = 1260. The resulting 
' suppression depends on the particle rapidity density dN/dy but weakly on the initial temperature. 

The jet energy loss treatment is also included in the calculation of high pr photons. Photons coming 
from primordial hard N-N scattering are the dominant contribution at RHIC for pr > 5 GeV, while 
£NJ ■ at the LHC, the range 8 < pr < 14 GeV is dominated by jet-photon conversion in the plasma. 

5h ! 

I. INTRODUCTION 

ON ' ....... „ 

' One of the most important issues that arises in the study of relativistic heavy ion collisions is that of the creation of 

a quark-gluon plasma (QGP). From lattice QCD calculations, there is strong evidence that at around a temperature 

C*") ' of 170 MeV, or equivalently an energy density of about 1 GeV/fm 3 p|, there should be a phase transition between 

an ordinary confined hadronic matter phase and a QGP phase, with an associated change in the relevant degrees of 

22 ' freedom. At the Relativistic Heavy Ion Collider (RHIC), measurements support the assertion that the initial energy 

density of the created system reaches up to 5GeV/fm . However, this does not necessarily imply the creation of 

■ a QGP. To make sure that a new state of matter has been formed, independent evidence is needed. One of the 

most striking measurements in support of the creation of hot and dense matter at RHIC is the discovery of high 

|/~} i pt suppression in central Au-Au collisions. This phenomenon is observed in single hadron spectra Q,0| and in the 

' disappearance of back-to-back correlations of high pt hadrons |j] . 
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There are a number of factors that can potentially influence the spectrum of high pt partons in heavy ion collisions 



O ■ compared to that in hadron-hadron collisions: 

(i) A difference between the parton distribution functions of a proton and a heavy nucleus. This can be both 
depletion (shadowing) or excess (anti-shadowing) depending on the value of momentum fraction x. This also 
includes the effect of gluon saturation. 



(ii) The initial state multiple scattering effect. This is the well known Cronin effect caused by multiple soft scatterings 
a parton may suffer before it makes a hard collision. This can, of course, only be significant in interactions 
involving nuclei. 

(iii) Final state energy loss. This is due to the interaction between the produced hard parton and the hot and dense 
environment. 

Our kinematical range of application in this work is mainly mid-rapidity. There, initial state effects cannot explain 
high pt suppression, as otherwise such suppression should also be observed in d-Au collisions, which is not the 
case. High px suppression has to arise from a final state effect: jet energy loss[5|. On top of experimental evidence 
for the suppression cited above, it has also been proposed [|| that azimuthal anisotropy at high px could also 
be explained by jet energy loss. Induced gluon bremsstrahlung, rather than elastic parton scattering, has been 
identified to be the dominant mechanism for jet energy loss [a 0. In the thermal medium, a coherence effect, 
the Landau-Pomeranchuk-Migdal (LPM) ^(j effect, controls the strength of the bremsstrahlung emission. Several 
models of jet quenching through gluon bremsstrahlung have been elaborated: Baier-Dokshitzer-Mueller-Peigne-Schiff 
(BDMPS) [UGl, Gyulassy-Levai-Vitev (GLV) 0, Kovner- Wiedemann (KW) and Zakharov There have 
also been phenomenological studies where different energy loss mechanisms were added onto Monte Carlo jet models 
fl6l fl7i flSl| (see al so Jl9l for a recent review). In this paper, we will use the formalism developed by Arnold, Moore, 
and Yaffe (AMY) 20], which correctly treats the LPM effect (up to 0(g s ) corrections), to get the parton energy loss. 
That model will be convolved with a (1 + 1)D expanding QGP to obtain quantitative results on jet quenching. Then, 
the electromagnetic signature of jet-plasma interactions will be investigated. 




This paper is organized as follows. In Sec. II we give a qualitative description of gluon radiation and briefly review 
the LPM effect and the conditions for the LPM regime. In Sec. Ill, the AMY formalism is presented and the main 
distinctions between the different treatments are highlighted. In the following sections, the complete space-time 
evolution of the jet is calculated within a hydrodynamical model, and high pr spectra are obtained for pions (Sec. IV) 
and real photons (Sec. V), for both RHIC and Large Hadron Collider (LHC) energies. Finally, Sec. VI contains a 
summary and conclusions. 



II. QUALITATIVE ARGUMENTS 

The high px data from various RHIC experiments (see Ref. [2lJ for a recent review) can be characterized by the 
following quantity: 

RAA = I AT \(JKT I A M V ~ C fOT PT ~ 2 GeV ' 1 

( A'coll) (dNp P /dyd z p T ) 

which represents the ratio of the number of events per unit rapidity and transverse momentum, scaled to the number 
expected based on proton-proton rates and the number of nuclear collisions occurring in the heavy ion collision. As 
indicated, this ratio is found experimentally to be approximately constant, with the constant C ranging from about 
0.2 to 0.6 depending on the centrality class and the rapidity y. The pp spectrum at mid-rapidity can be conveniently 
parametrized by |22j 

dN pp a ( Po 



= A l^_\ , (2 ) 
dyd z pr \Pq+PtJ 

where A is a normalization constant, po — 1.2 GeV and n — 10. In a previous paper two of us showed that the 
leading order thermal QCD result, shown in Fig. ^ can qualitatively explain the fact that C is roughly constant over 
a finite kinematical window. In this section, we use the simple gluon radiation spectrum used by Baier et. at. [Tl) to 
illustrate the main physical features of this calculation. 

As previously stated, the dominant energy loss mechanism in medium is the radiation of gluons. We will express the 
rate of radiation of gluons of energy w, per unit time and u> interval, as T = dl/dwdt. To understand the relationship 
between the behavior of Raa and T, we note that, given T, the parton distribution function P(p,t) (which gives the 
probability that there will be a parton of energy p at time t) should satisfy a rate equation, 



dPfat) 
dt 



dojP(p + uj,t)r(p + uj,uj,t)- P{p,t) / dujT(p,oj,t), (3) 
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provided that conversion between different species of partons is unimportant. If the rate T(p, u>, t) has a weak depen- 
dence on p (an assumption we make here for illustration but will not use in our work in subsequent sections), then 
an approximate solution can be written as 



P(p,i) w / deD(e,p,t)P (p + e) 



(4) 



where 



D(e,p,t) = exp 



dui J(p, u>, t) 



oo ^ 

71 = 



n 



i=l 



with 



J{p,uj,t) = / dt'T(p,uj,t') . 



(5) 



(0) 



A similar expression is written down in Ref. 12]. The main difference between the current expression and the one 
in |12| is that the range of integration in the above equations is from — oo to oo: The negative uj corresponds to the 
absorption of thermal partons and the range uj > p corresponds to annihilation of a parton with a thermal anti-parton. 
Following we let 



Po(p + e) 
Po(p) 



exp(-nefp) , 



(7) 



which allows us to write the nuclear modification factor as 



Raa{p) = J deD(p,e,t) 



Po(p + e) 



exp 



Mp) 

J de D(p,e,t) exp(— ne/p) 

/oo 
-oo 

For this to be approximately constant, or at least have a weak dependence on p, the exponent 

/OO pOO 
dujJ(p,uj,t) + / dwJ(p + W , W ,i)e-" /(,,/,l) 
-oo J — oo 



(8) 



(9) 



must have weak dependence on p, at most equaling /31np with /3 < 1. 

Even without knowing the explicit form of the J function, some insight can be gained from this expression. First, 
for the second integral to be well defined in the u> < range, the Boltzmann factor corresponding to the thermal 
partons must compensate the e~"^ p /"^ factor. Secondly, because of the e^^l ( -Pl n ) factor, in general the loss term 
dominates and hence Raa(p) < 1- Third, as p — * oo, the compensation between the two terms in A becomes complete 
and Raa — * 1- 

To be more specific, the precise form of the J function needs to be known. The energy loss of an energetic parton in 
a dense medium mainly proceeds via bremsstrahlung of soft gluons. The form of the energy loss (the J function above) 
depends crucially on whether bremsstrahlung off many scatterers are coherent or incoherent. This is determined by 
the relative size of the mean free path A = 1 /na and the coherence length [2J| , 

(10) 



where fj, is the typical size of the soft momentum exchange and to is the energy of the emitted gluon. When Z co h <C A, 
then one is in the extreme Bethe-Heitler regime where the energy loss per unit length is proportional to the incoming 
energy. On the other hand, when Z co h 3> A, then one is in the extreme Landau-Pomeranchuk-Migdal (LPM) regime 
where the energy loss per unit length is proportional to the square-root of the incoming energy. 
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The typical mean free path for soft scattering in a hot medium can be estimated as follows. The density of scatterers 
according to the thermal distribution is n ~ T 3 . The typical soft (p e xch ~ gT) scattering cross-section is given by 
cr - g 2 IT 2 . Hence 

\~l/g 2 T. (11) 
According to the above estimate of the coherence length, the LPM effect becomes relevant when ^ co h > A or 

- 2 > A, (12) 
M 

with u ~ gT and A ~ l/g 2 T. This yields the condition 

u>T. (13) 

Defining -Elpm = A/x 2 , one can then conclude that for the emitted gluon energy in the range of < lo < -Elpm, energy 
loss is governed by the Bethe-Heitler limit 

dl ct s ^ 1 (14) 
diodz it c A ' 

and for lo > -Elpm, the energy loss is governed by the LPM limit |24j . 

(15) 

For lo > Blpm(-^/A) 2 , the coherence length exceeds the length of the medium. In this case, effectively only a single 
scattering can occur and the radiation spectrum goes back to 

lo-^—~—N- (16) 

In this paper, we will only consider the case 

< ElPM < Eo < -Ef ac t (17) 

assuming that L/X^> 1. Here Eq denotes the original energy of the parton and we defined -Ef ac t = -Elpm^/A) 2 . 

In order to use Eq. ©, in addition to the above forms of gluon spectrum for lo > 0, we also need to know the form 
of J(p,ui,t) for uj < 0. For simplicity, we take the spectrum for lo < to be the same as the Bcthc-Hcitler spectrum 
and multiply it by the Boltzmann factor to take into account that absorption of high energy gluons from the thermal 
medium is Boltzmann suppressed: 

dl a,iV„ 1 i i irp , 
= _^£ e-H/r. 18 

diodz 7TA \uj\ 

For illustration, we set -Elpm = 1 GeV, A = 1 fm and T = 400 MeV and evolve the system until L = ct = 4 fm. 
We have also set the upper limit of the momentum integral to p. The above simple gluon spectrum then yields the 
Raa(p) curves shown in Fig[3 The solid curves in the figures are calculated including the absorption contribution 
and the broken curves are calculated without the absorption part. Also shown are the dot-dashed curves which are 
calculated with only the lo > £?lpm part of the gluon spectrum omitting the Bethe-Heitler part. In Fig|3J we show a 
momentum range far exceeding the factorization energy to compare with the calculation in Ref. [23| . 

In Ref |28j | , two of us have calculated the same ratio of the final and the initial spectrum using the resummation 
method developed by AMY. The shape and the trend of the curve calculated with the thermal absorption in Figs [21 
and [3] closely match up with what we have previously calculated in Ref |23j. The experimentally obtained Raa from 
STAR and PHENIX from RHIC experiments also roughly has the same shape and trend EsT]. 

There are two major points we would like to make with this simple calculation. One is the importance of the 
absorption of the thermal gluons. This is also closely related to the appearance of a momentum scale p = nT. Near 
this scale, absorption of thermal gluons is very efficient in changing the shape of the final spectrum. Comparing the 
two curves in Fig|21and also in Fig|31 it is clear that the roughly flat shape within the limited momentum range shown 
in Figgis due to the absorption of thermal gluons. Furthermore, the efficient absorption near p — nT causes Raa to 
have a larger value near this point. Without the absorption, Raa is a monotonically increasing function of p. In the 
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FIG. 2: Approximate -Raa(p) calculated with Eq. @ following the approach of BDMPS up to the maximum momentum 
of -Efact = 16 GeV. The solid line includes the contribution from absorption of thermal gluons whereas the broken line only 
includes the emission. The dashed-dotted line is calculated with only the LPM part of the gluon spectrum. The divergence at 
p = 10T is an artifact of the approximation. 
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FIG. 3: The same as Fig|5]but showing a larger energy range. 

experimentally measured Raa, one does observe such a change in slope near px — 3 GeV. This is usually attributed 
to the Cronin effect. However, it could very well be due to the absorption of thermal gluons. 

The importance of gluon absorption will arise whenever the slope of the spectrum, dhiP/dp, is comparable with 
the thermal value 1/T. When these are equal, detailed balance ensures that the more numerous lower energy particles 
keep re-populating the higher energy ones. The reason the spectrum in Fig. prises to 1 at 10T is because we took the 
initial particle distribution to be a strict power law, ~ p~ n , which becomes steeper than the thermal spectrum below 
p = nT. More realistic initial distributions, as in Eq. will not display this behavior; nevertheless, absorption is 
important in the lower energy region. 

Another point to make is the importance of the Bethe-Heitler part of the radiation spectrum. The effect of neglecting 
the Bethe-Heitler part of the spectrum becomes justifiable only around p ~ 100T. Therefore, calculations with only 
the LPM part of the spectrum appear unrealistic for RHIC energies. 
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FIG. 4: A typical bremsstrahlung diagram that needs to be resummed. 




FIG. 5: A typical gluon self-energy diagram that needs to be resummed. 



III. THE FORMALISM 



In a seminal series of papers, Baier, Dokshitzer, Mueller, Peigne and Schiff identified the diagrams that contribute 
to the leading order energy loss in a hot QGP medium. A typical diagram is shown in FigQ] To compute the energy 
loss rate, one must deal with the following issues: 

1. How are these diagrams to be resummed? 

2. How is the resulting equation to be solved? 

3. How does one describe the time evolution of the energy loss? 

We now summarize the approach taken by Arnold, Moore and Yaffe (AMY). The ingredients discussed in this section 
can be found in Ref. |23j. We shall also briefly compare with other approaches. 

In the diagram shown in Fig. 0] a single wiggly line connected to a single (g> is not a gluon propagator. Rather, it 
just denotes that a quark has interacted with the soft background field at a certain time ti or Si. To get the gluon 
energy spectrum at the final time, we need to sum all the diagrams of the type shown in Fig. 0] then square and 
average over the soft field configurations. This procedure leads to the conclusion that the scattering rate can be 
obtained by calculating the imaginary part of the gluon self-energy in the soft thermal background. A typical gluon 
self-energy in such a situation is shown in Fig. [5J 

The result of resummation is summarized in the following formulas. The transition rates for various emission 
processes are given by 



dT(p,k) 



Cs9 2 s 1 1 



l+(l-x) 2 
x 3 (l— x)' 2 




dkdt 



x 




X 



G 



-^2h-ReF(h,p,fc). (19) 

Here C s is the quadratic Casimir relevant for the process (in QCD, 4/3 for processes involving a quark and 3 for the 
pure glue process), and x = k/p is the momentum fraction in the gluon (or the quark, for the case g — > qq). The 
factors l/(l±e- k / T ) are Bose stimulation or Pauli blocking factors for the final states, with — for bosons and + for 
fermions. h = p x k determines how non-collinear the final state is; we treat it as parametrically 0(gT 2 ) and therefore 
small compared to p • k. Therefore it can be taken as a two-dimensional vector in the transverse space. F(h,p, k) is 
the solution of the following integral equation: 

2h=tSE(h,p,k)F(h)+g 2 J T^kC(q ± ) x {(C s - C A /2)[F(h) - F(h-fccu)] 

+ (C A /2)[F(h)-F(h +J3 q ± )] 
+(C A /2)[F(h) - F(h-(p-fc)q ± )]} • (20) 

Here C(qj_) is the differential rate to exchange transverse (to the parton) momentum q_L. In a hot thermal medium, 
its value at leading order in a s is |26| 



qiCqi+mS,)' "* u ~ 6 



2 7Z2 Z^T > mi ) = ^_(2iV c +JV f ). (21) 



The energy difference between the final and the initial states is given by 

SE (h v k) h2 + < + m P- k m ? (22) 

dE[h ' P > k) 2pk(p-k) +1 2k + 2fr=jfe) ~~ 2^ ' (22) 

where m 2 arc the medium induced thermal masses, equal to m 2 D /2 for a gluon and C/<? 2 T 2 /4 — g 2 T 2 /3 for a quark. 
For the case of g — > qq, (C s — Ca/2) should appear as the prefactor on the term containing F(h — pqj_) rather than 
F(h-fcq ± ). 

Next, we use these expressions to evolve the hard gluon distribution P g (p,t = 0) and the hard quark plus antiquark 
distribution P J q{p 1 t — 0) with time, as they traverse the medium. The joint equations for Pjq and P g are 



where the k integrals run from — oo to oo. The integration range with k < represents absorption of thermal gluons 
from the QGP; the range with k > p represents annihilation against an antiquark from the QGP, of energy (k—p). 
In writing Eq. (|2*3^l . we used dT^ g (p,k) — dT^ g (p,p—k) and similarly for g — > qq; the 6 function in the loss term 
for g —* gg prevents double counting of final states. Since brcmsstrahlung energy loss involves only small 0(g s T/p) 
changes to the directions of particles, these equations can be used for the momentum distributions in any particular 
direction. For a single initial hard particle, they can be viewed as Fokkcr-Planck equations for the evolution of the 
probability distribution of the particle energy and of accompanying gluons. These expressions depend at several points 
on g 2 or a s . When evaluating them numerically, we have used a s = 0.3. 

We now briefly compare our approach with others used in the literature. As this paper is not intended to be a 
review, this cannot and will not be a complete comparison with all other jet energy loss calculations. 

As was mentioned at the beginning of this section, most jet quenching calculations start with the problem of 
resumming all diagrams of the form shown in Fig0] In a series of papers, U. Wiedemann et al. have shown that 
BDMPS, Zakharov, GLV and also the eikonal approach taken by Wiedemann and Kovner are all related to each 
other. In these approaches, the problem of resumming the diagrams is solved in position space assuming static 
scatterers, although the methods of solving the problem differ. Since a thermal medium consists of dynamic, moving 
scatterers, temperature was introduced only as the controlling variable for the mean free path in these approaches. 
The evolution of the initial parton distribution is then achieved through Eq. with the important restrictions that 
only the emission of gluons was considered with (Wiedemann) or without (BDMPS) restriction on the kinematic 
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upper limit of the emitted energy. There is also an approximation that ignores the emitted energy in the gain rate 
(T(p + u,u,i) -> T(p,cj,t) in Eq. ©). 

Our approach differs with the above ones in three major ways. The biggest difference is that our calculation is 
completely thermal and hence the scatterers are all dynamic. In our calculation, temperature enters through the 
thermal phase space of the initial and the final particles and there is no assumption of the form of the elementary 
cross-section. All are calculated completely within the framework of hard thermal loop resummed leading order 
thermal QCD. Hence, gain or loss due to the absorption of thermal partons (u) < 0) as well as the loss process 
of pair annihilations with the thermal partons (u> > p) are fully included in our calculation. These are missing in 
the approaches mentioned above. In Ref. jS^, thermal absorption and stimulated emission were introduced in the 
framework of GLV. This was done only up to the first order in the opacity expansion without including the annihilation 
process with thermal partons. Second, instead of approximately solving for the transition rates and then using Eq. JIJ, 
we explicitly solve for both the transition rates and the coupled rate equations of hard quarks and gluons. There 
are also no approximations about whether the Bethe-Heitler or LPM regime is relevant; the transition between these 
extremes is handled correctly. 

One limitation of the current approach as compared to the previous ones is that the transition rates are calculated 
in momentum space assuming the thermodynamic limit. That is, our approach assumes that the high energy parton 
experiences a uniform medium on the time scale of the formation time of the emitted radiation. Hence, our approach 
is limited to momenta less than the factorization energy. With fj, w 0.5 GeV, L « 5 fm and Awl GeV, this limits the 
momenta to the region p < 30 GeV. This is not a big problem at RHIC energies, where y/s is only 200 GeV, and it 
also covers the px acceptance of the ALICE detector at the LHC. 



IV. PION PRODUCTION 



The goal of this section is to use the formalism explained in the previous section to calculate the neutral pion 
spectrum in heavy ion collisions. Our approach to this problem relies on the fact that for hard spectra, the AA 
collision can be regarded as a collection of binary collisions. In this way of formulating the problem, the AA spectrum 
is given by the convolution of the elementary pp spectrum with geometrical and in-medium factors. 

Up to suppressed corrections, the 7r° cross-section in an N-N collision factorizes and is given by |28j 

E n — 3 = V" dx a dx b g(x a ,Q)g(x b ,Q)K jet — D n o/ c (z, Q ) (24) 

a Pn J ctt 7TZ 

a.b.c.a 

where g(x, Q) is the parton distribution function (PDF) in a nucleon, D n o i c (z, Q ) is the pion fragmentation function, 
da + dt — — is the parton-parton cross-section at leading order and the factor Kj et accounts for higher order effects 
(where "jet" here means a fast parton having pi?* 3> 1 GeV). According to j2^, Kj et is almost p^f* independent at 
RHIC energies. We will use Kj et ~ 1.7 for RHIC and 1.6 for the LHC, based on their results. 

For all our calculations, we set the factorization scale (Q) and the fragmentation scale (Q ) equal to pr- We take 
the CTEQ5 parton distribution function |30| and the it fragmentation function extracted from e + e~ collisions |3l| . 

Fig. H3 shows our calculation for the spectrum of high px neutral pions in pp collisions at RHIC, compared to the 
PHENIX result (3^1 ■ O ne can readily see that our calculation reproduces the data in the region where jet fragmentation 
is expected to be the dominant mechanism of particle production (pt > 5 GeV/c) |33j. 

To obtain the high pt 7r° cross-section in AA collisions, we must modify the pp calculation in two ways. First, the 
PDF of a nucleus differs from that of a proton; 

gA(x a ,Q) = g(x a , Q)RA(x a ,Q) , (25) 

where the nuclear modification of the structure function Ra takes into account shadowing and anti-shadowing. In 
this paper, we use RA(x a , Q) as parametrized by Eskola et al. |34j |. 

Second, we must account for the energy loss of the parton between its production in the hard initial scattering 
event and its hadronization. We will assume that a jet fragments only outside the medium, as may be justified by 
estimating the formation time of a pion with the typical observed energy. In Ref.[35|], the formation time of an emitted 
pion with energy E v is estimated to be ~ R^^-, where R^ ~ 1 fm is the size of the pion. For a 10 GeV pion, this 
gives Tf ~ 35-70 fm/c. This is much longer than the path of a jet in the hot matter created by two colliding nuclei. 
Therefore, one should consider production, energy loss in medium, and fragmentation to occur sequentially. 

We will assume that the fragmentation at the edge of the QGP involves the usual vacuum fragmentation function. 
The medium effect is then to reduce the parton energy by an amount determined by the Fokker-Planck equation 
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FIG. 6: (Color online) Neutral pion spectra in pp collisions at RHIC. The data points are from PHENIX and the solid line is 
the calculated result from jet fragmentation. 



previously presented, Eq. (|23|l . This is most conveniently written by denning a new medium-inclusive effective 
fragmentation function, 



Ar°, c (z,Q;r,n) = / dpf—(P qq/c (pf : p t )D^a /q (z l ,Q) + P g/c (pf,pi)D 7I o /g (z',Q)) , (26) 



z 
z 

where z — pr/Pi and z' — px/pf- Pqq/ c {PflPi) and Pg/JjPf'iPi) represent the solution to Eq. I)23|l . which is the 
probability to get a given parton with final momentum pf when the initial condition is a particle of type c and 
momentum pt. 

The quantities P q q/ C and P g i c in Eq. I|26|) depend implicitly on the path length the initial parton must travel through 
the medium and the temperature profile along that path. This is not the same for all jets, because it depends on 
the location where the jet is produced and on the direction the jet propagates. Therefore, one must convolve this 
expression over all transverse positions r_L and directions n. Since the number of jets at r_L is proportional to the 
number of binary collisions, the probability is proportional to the the product of the thickness functions of the colliding 
nuclei at rj_. For central collisions where the impact parameter b ss 0, we get 

V(t ± )<xT a (t ± )T b (t±). (27) 

For a hard sphere which we use for simplicity, this probability is 

2 



PCrxH^l-^j^-rO, (28) 

which is normalized to yield J d 2 r± V(r±) = 1. Since the direction of the jet is fixed by the pion direction (n = rjpj), 
the final in-medium modified fragmentation function is 



D n0/c (z, Q)=l d 2 r ± P(r ± )D n0/c (z, Q, rj_, n) . (29) 
The AA spectrum is now given by 

dN AA (Neon) f , , ( n \ l n\v da a + b ^+ d D„o /c (z, Q) 

-7—^ — = > dx a dx b g A (x a ,Q)g A (x b ,Q)K jet — , (30) 

dyd A p T a in r^.J dt nz 

a,b,c,d 
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FIG. 7: The evolution of V(AE\E C ) as a function of time (passage length). 



where (N co \\) is the average number of binary collisions and (Tin is the inelastic nucleon-nucleon cross-section. 

Since Eq. (|30[1 is expressed in terms of probability distributions, it is straightforward to evaluate it using the Monte- 
Carlo method. One complication in doing so is that, while the parton traverses the medium, the medium also evolves. 
Therefore at each time step of solving Eq. (|23[) . the temperature must be adjusted to the local environment. We 
assume that the medium expands only in the longitudinal direction, based on the following reasoning. The low pt 
neutral pion spectrum at RHIC is well reproduced by a hydrodynamical model incorporating transverse expansion, 
while the model fails for p? > 3 GeV, suggesting that high-p-^ pions mainly come from jet fragmentation 36]. The 
transverse expansion will have two effects on jet energy loss. First, the expanding geometry will increase the duration 
of parton propagation. However, the same expansion will make for a falling parton density along the path. Those two 
effects partially compensate each other and the energy loss is just about the same as in the case without transverse 
expansion [37| . In 1-D Bjorken expansion [ssj . the temperature evolves as 



T = Ti 



(?) 



1/3 



(31) 



In the original Bjorken model, the transverse density is assumed to be constant. Since a nucleus does have a 
transverse density profile, it is more realistic to assign the initial temperatures in the transverse direction according 
to the local density so that |3^, |4£j 



T(r,n) = T, 



2 1- 



Ri 



-,1/4 



Putting Eqs.JSU an d l|32(l together, we get the temperature evolution of a QGP expanding in 1-D as 



T(r,r) = T i (^) 



1/3 



2 1 



1/4 



(32) 



(33) 



The jet evolves in the QGP medium until it reaches the surface or until the temperature reaches the transition 
temperature T c . In our calculation, we assume a first-order phase transition and use 



Iqgp 



1 



1 



(34) 



as the fraction of the QGP phase in the mixed phase |3^| . Here = j^- is the ratio of the degrees of freedom in the 
two phases and Tf is the time when the temperature reaches T c . The evolution equation is then scaled accordingly 
for t > Tf. We take the critical temperature T c to be 160 MeV. 

The result of our calculation for RHIC energy is summarized in Fig. [S] together with the PHENIX data Q ■ To 
cover the uncertainties in the initial conditions, we consider two different sets, one at a relatively high temperature 
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FIG. 8: (Color online) Nuclear modification factor for pions at RHIC. Data points are from PHENIX 0. The solid lines show 
the full calculation of the spatial distribution of jets in the plane z — for two initial temperatures. The dotted lines assume 
that all jets are created at the center and the dashed lines assume the same approximation but with a reduced energy loss. 



of Tj = 447 MeV and a relatively short initial time of Tj = 0.147 fm/c taken from Refs. H^,^3|) an d one at a lower 
temperature Tj = 370 MeV and somewhat later time for the hydrodynamic evolution of Tj = 0.26 fm/c taken from 
Ref. |4l|. Those two sets correspond to dN/dy\ y —o ~ 1260, estimated for central collisions in Ref. j4J|. In deriving 
these results, we have used a s = 0.3. If the value is larger then energy loss will be greater, so Raa will be smaller; if 
it is smaller, then Raa increases. Besides this dependence on a s , our results rely on no free parameters. 

The solid lines in the figure are full calculations with the initial spatial distribution V(r±) given in Eq. I)28[l. For 
comparison, we also display two more sets of calculations. The dotted line is calculated with the jets originating only 
from the center of the disk (rj_ = 0). Comparing the dotted line with the solid line, it is apparent that the absolute 
magnitude of the Raa depends very much on the density profile of the nucleus. The dashed line is calculated again 
with the jets from the center but with the energy loss rates reduced by a factor of 0.64. One may say that the average 
path of a jet has the length of about 0.64 x Ra- 

In Fig. [5J one can see that both Tj = 447 MeV and Tj = 370 MeV describe the real data reasonably well. This is 
somewhat surprising. Since the density of thermal particle is proportional to T 3 , the density at 447 MeV is about 
1.8 times the density at 370 MeV. Yet the energy loss does not reflect such a big difference. The reason is because 
the energy loss depends mostly on the duration of evolution. In a Bjorken expansion, the initial time Tj and the 
temperature are related by 

3 7T 2 1 dN 

Hence once dN/dy is fixed, the time evolution of the temperature follows a common curve regardless of the initial 
temperature (c.f. Eq. I|33|) ). The only difference between the higher and the lower temperature cases is that the time 
evolution starts earlier for the higher temperature case. From the moment r passes Tj for the lower temperature, the 
evolution of the two systems is identical. Therefore, if energy loss at the beginning of the evolution is small compared 
to the later time energy loss, the amount of energy loss depends mostly on the duration of the QGP phase At — Tt—Ti. 
Since r/ S> Tj, the duration of the QGP phase is approximately the same for high temperatures. We have verified 
that the energy loss between the times corresponding to Tj = 1000 MeV, T, = 447 MeV and Tj = 370 MeV are at most 
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FIG. 9: (Color online) Nuclear modification factor for pions at RHIC, taking all jets to originate at the center of the nucleus. 
Data points are from PHENIX |2|. The solid lines are for initial conditions which lead to a particle multiplicity of dN/dy = 1260 
while the dotted lines are for initial conditions leading to dN/dy = 680. For each set of lines, the initial temperature is, from 
bottom to top, 1000, 447 and 370 MeV. 



about 10%. 

Since the suppression is mainly controlled by the duration of the evolution, the suppression should be sensitive to 
the particle rapidity density dN/dy, which fixes the lifetime of the QGP. We can see in Fig. EH that it is indeed the 
case. For simplicity, we assume here that the jets are all created at the center of the system. We see that there 
is very little change in the suppression as the initial temperature varies from Tj = 370 MeV to Ti = 1000 MeV. 
However, the suppression shows a strong dependence on dN/dy; going from dN/dy = 680 to dN/dy — 1260 increases 
the suppression by a factor ~ 1.5. 

For the LHC, we use the initial temperature from Ref. 41], T t =845 MeV, giving r 4 = 0.088 fm/c for dN/dy ~ 
5625 H3, H3- At this energy, jets could in principle be produced with energies as high as y/s/2 = 2750 GeV. 
However, according to our discussion in Sec. [HI the contribution from a jet having more that twice the observed pion 
energy should be sharply cut off by the steeply falling initial function. Since our approach is limited to the observed 
Pt 30 GeV, considering original energy of only up to about 100 GeV should be enough for our purposes. To be on 
the conservative side, we cut off the maximum jet energy at 400 GeV. 

The calculated LHC nuclear modification factor is shown in Fig. EH with the full fragmentation function from 
Eq. (|29[l . We also show in the same figure the effect of not tracking the secondary partons (gluons emitted by quarks 
and quarks emitted by gluons). Comparison between the solid line and the dashed line shows that ignoring the 
secondaries can make an overall difference of 30%. This arises almost entirely from quark secondaries from gluon 
jets-that is, from hard gluons which split, due to plasma interactions, into qq pairs, which subsequently fragment. 
This is important at the LHC, but much less so at RHIC, because most hard jets at the LHC arise from gluons, 
whereas a larger fraction at RHIC are from quarks. The gluons emitted by quarks are mostly soft; they also lose 
energy quickly and fragment inefficiently. At RHIC energies the error from dropping secondaries is only 5%. 

Comparing the RHIC and LHC results, we see that the suppression due to jet energy loss at RHIC is about a factor 
of 3 for a 10 GeV pion, while the suppression at the LHC for a pion at the same energy is about a factor of 6. This 
difference arises because the medium at the LHC remains hot for much longer than at RHIC. The cooling is governed 
by the product T 3 r oc dN/dy. Hence, smaller dN/dy implies faster cooling. We have also calculated the impact of 
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assuming a first order phase transition. At RHIC, at the end of the QGP phase, a large number of jets has already 
left the medium, but those still inside will suffer additional suppression during the mixed phase, such that assuming 
a cross-over between the QGP and the hadron gas phase will increase Raa by ~ 20%. At the LHC, at the critical 
point, a much larger fraction of jets are out of the medium, such that Raa is quite insensitive to assumptions related 
to the transition. Our initial conditions correspond to an initial time smaller than the common value of Tj= 0.6 fm/c 
used in hydrodynamic calculations j^. We have verified that going from t; = .26 fm/c (Tj = 370 MeV) to 0.6 fm/c 
(Ti=280 MeV) increases Raa by less than 20%. In this section, the AMY formalism has been applied to successfully 
reproduce the tt° spectra at RHIC. We will see in the following section, how the same formalism can be applied to 
determine the production of real photons. 
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FIG. 10: (Color online) Nuclear modification factor of pions at the LHC. The solid line includes pions coming from 
bremsstrahlung secondaries emitted in the thermal medium; the dashed line does not. For both lines, the full spatial dis- 
tribution from Eq. I28H has been used. 



V. PHOTON PRODUCTION 



The hard photons produced in nucleon-nucleon collisions can be divided into three categories: direct photons, 
fragmentation photons, and background photons. Direct photons are those produced by Compton scattering and 
annihilation of two incoming partons. Fragmentation photons are those produced by bremsstrahlung emitted from 
final state partons. Background photons are those produced by the decay of hadrons subsequent to the collision, 
primarily from 7r° — * 77 decay. The 'prompt photons' are those coming from direct production and the fragmentation 
process. The expression for prompt photon production is 

E ^ = E J dx a dx b9 (x a ,Q)g(x b ,Q) {lf 7 (p T )- ^ ^f^S (x h - X *$**^J) 

+K brem (p T ) — D~/ c (z,Q) > . (36) 

at ttz J 

K 1 and K brem are correction factors to take into account NLO effects; we evaluate them using the numerical program 
from Aurenche et al. 0, obtaining K 7 {10 GeV)~ 1.5 for RHIC and LHC and K brem (10 GeV)~ 1.8 at RHIC and 1.4 
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FIG. 11: (Color online) Prompt photons produced in pp collision at RHIC. Data points are from PHENIX. The solid line is 
calculated with Eq. 136H . 



at LHC. All scales (renormalization, factorization and fragmentation) have been set equal to the photon transverse 
momentum px- The photon fragmentation function D^/ c is extracted without medium effects in e~ +e + collisions |45l |. 
The validity of this expression for pp collisions at y/s = 200 GeV is shown in Fig^Jwith data from PHENIX (Klaus 
Reygers, Hard Probes 2004). It appears clear that the baseline mechanism of high px photon production in nucleon- 
nucleon collisions is under quantitative control. 

In AA collisions, there is an additional source of high px photons: the medium contribution. This contribution 
includes the direct conversion of a high energy parton to a high energy photon by annihilation with a thermal parton, 
in-medium bremsstrahhmg from a jet, and thermal production of photons. 



A. Jet-Photon Conversion 



In Ref. |39j, the conversion of a leading parton to a photon in the plasma was found to be an important process. This 
happens when a jet crossing the hot medium undergoes an annihilation (q + q — > g + 7) or a Compton (g + q — ► q + 7) 
process with a thermal parton. The related production rate of photons is given by |39l |46| 



dR sr^ I ef\ 2 T 2 aa s . „ ... „ , ... 

E T -^If^ + fM)] 



dycPpx ~* V e / 8tt 



2 111 I ^— I — Cann — Ceo 



(37) 



where T is the temperature, C ann =1.916, and Cc O m=0-416. It has been shown in Ref. (50l l51j. that the infrared 
dependence on the quark mass should be replaced by a dependence on the in-medium thermal quark mass, m 2 — 

d 3 xd 3 p 



2m 2 h = Aira s T 2 /3. The phase-space distribution function of the incoming particles is denned as 

9i I ^ff(x, P ) = N l , (38) 



where Ni is the number of particles i and gi is the spin-color degeneracy. The phase-space distribution function for 
an incoming jet, assuming a Bjorken -q — y correlation [52|. is 

fjet{z,P,to) = -7-75 — o{r)-y) 

g q rp T dyd 2 p T 

(27rfr(v ± )t p T dN^ 

= m 9 ^2l~n — o(z-w^ ), (39) 

g q \Jt 2 - z 2 E 2 dyd 2 p T 
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where r] is the space-time rapidity and to is the formation time of the jets. As before, the jets are taken to be massless 
and we suppose that energy loss of jets in the plasma does not change their direction. With this latter approximation, 
fj e t can be factorized into a position space and a momentum space part: 



1 dNi et (E,t) 
E 2 



dE 



f f 1 dW*(E,t) 



(40) 



where t — t — to is the propagation time of the jet. In the high energy limit, the cross-sections a q q^ gl and &qg^ q ~t 
are dominated by direct exchange between the quark and the photons. Since we are interested in photons produced 
at mid-rapidity (y — 0), we only need to consider quark and anti-quark jets produced at mid-rapidity. This gives 



fjet(~&,T?,t)\ 



(2Tr) 3 V(\wl\)t 1 djyw 
lv= ° W*o ~ z 2 PTdyd 2 pT 
where w$ is, in the plane z = 0, the initial radial position of the jet, 



(pT,t)6(z) , 



= \-5t-t-, 



P 



(41) 



(r cos 4> - t) 2 + r 2 sin 2 (f> for t ~ , (42) 

and cj> is the angle in the plane z — between the direction of the photon and the position where this photon has been 
produced. The AN 
is Iqq(E, t), where 

dN qq dN™ 
P qq (E,t) = - PEr (E,t)<xp T ^ : ^ z -(p T ,t). (43) 



produced. The AMY formalism is introduced here to calculate the evolution of J^p pT (pti £)• The output of Eq. I|23|) 



dE 



dyd 2 p T 



The initial distributions Iqq(E,to) and P g (E,to) are fixed by the initial jet distribution at mid-rapidity as 
parametrized in Ref. 39, 40] . The total photon spectrum is given by a full space-time integration: 



dNi 
dyd 2 pj 



= / drr I rdr I d(p I dq 



dR 



dyd 2 p T 



(E 7 = pr cosh(y - r))) 



dt / rdr 



dz 



dR 



dyd 2 p T 



(E 1 = pTCOsh(yo)) 



(44) 



The production rate is calculated in the local frame where the temperature is defined where the photon rapidity 
becomes yo = y — ?/• 

High-pr photons are emitted preferentially early during the QGP phase, when the temperature is at its highest 
point. Indeed, explicit hydrodynamic calculations show that the nuclear space-time geometry smoothly evolve from 
1-D to 3-D |4g • By the time the system reaches the temperature corresponding to the mixed phase in a first-order 
phase transition, the system is still very much 1-D |47|. For such a geometry, specific calculations suggest that 
the flow effect on photons and dileptons from the QGP is not large at RHIC and LHC for px > 2 GeV. Assuming 
again a 1-D expansion, we get the production of photons from jet- medium interactions from Eqs. (|37|l . (|41(l and (144ft . 
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where the temperature T evolves according to Eq. (|33fl . The 4> integration can be done: 
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FIG. 12: (Color online) Direct production of photons by jets in the plasma for Au-Au at RHIC and Pb-Pb at the LHC. For the 
solid lines, jet energy loss is included; for the dashed lines it is neglected. The initial temperature is = 370 MeV at RHIC 
and T;=845 MeV at the LHC. 



where 



M 



arccos 



t 2 - R\ 
2tr 



(47) 



Then, the final expression for the jet-photon production becomes 
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As before, we assume a first order phase transition beginning at the time r/ and ending at th = fdJf |36 
Tf we scale the production rate by /qgp (Eq. (|34|) ). such that 
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(48) 



After 



(49) 



The first term and second term correspond respectively to photons produced during the pure QGP and the mixed 
phase. Our results for RHIC and LHC, with and without energy loss are shown in Fig.^| Here also, a s = 0.3. We 
see that higher energy photons are more sensitive to jet energy loss: photons at 4 GeV are suppressed by a factor 
1.3 at RHIC, while 15 GeV photons are suppressed by a factor 1.6 due to jet energy loss. We find approximately 
the same suppression for the LHC. This suppression is much smaller that the one observed from Raa in the previous 
section. This is because the photons can be produced at any point in the hard parton's propagation through the 
medium, while the jet energy loss depends on the final parton energy. Hence, some of the photon rate arises from 
before, rather than after, the jet has lost much energy. 



B. Bremsstrahlung Photons 



Hard partons in the medium can also produce photons by bremsstrahlung when they scatter in the medium. The 

dkdt follows the same expression as JkSt 



photon bremsstrahlung rate dT jujf'^ follows the same expression as dT dk J f p ' k ' > in Eq. i|19|) . but with C s — > % 



Ca — > 0, 772 -y = and — * e . The photon bremsstrahlung distribution is given by 



dt J qq ^ ' dkdt 
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It is assumed here that the photon production rate is so low that the quark plus anti-quark distribution P^q will be 
unchanged by the photon emission. The photon distribution is finally convolved with the initial spatial distribution 
of jets to get the final spectrum of bremsstrahlung photons, 



dN. 



jet—br 



dN iPt - 



jet—br 



- ( d 2 r ± V(r ± )P 7 ( PT ,d), (51) 

>T J 



dyd 2 p T y =o Pt^(u = 0) dp T p T 

where rj_ is the position where the jet has been created, and d = d(r±) is the distance crossed by the jet in the plasma. 
The factor Cl(y = 0) corresponds to a d(j> and a dy integration around the plane y = 0. This factor can be absorbed 
in the definition of the initial distribution Bjq(pr,t = 0). 

Numerically, bremsstrahlung photons turn out to be subdominant to jet-photon conversion. This is because, while 
the rate at which such photons are produced is larger, they typically carry only a fraction of the jet's energy, while 
jet-photon conversion predominantly produces photons with the complete energy of the jet (hard parton). When 
folded against a steeply falling spectrum of jets, the process which produces the highest energy photons will dominate 
the final spectrum. 

Recently, Zakharov has also considered bremsstrahlung emission of photons from jets |4flj . His work accounts for 
finite-size effects in the high energy limit; in this respect it is more complete than our work here. However, our energy 
range of interest for this study, as discussed at the end of section ITTT1 is below the factorization scale, so finite size 
effects are small for the energies we consider. 



C. Thermal Photons 



The thermal-thermal contribution comes from the photons produced by two scattering thermal particles. The 
Compton and annihilation rate has been calculated in the previous literature [50l l5lj . and we use those results 
here. There are also leading order bremsstrahlung and inelastic pair annihilation contributions |53|: we use the 
parameterization for those rates presented there. 



D. Non-thermal contributions 



The expression for prompt photons produced in AA collisions is 



dNj-projnpt _ ^ da (N co u) 
dyd 2 p T 7 d 3 p 7 a m 

where we take the values (iV co u) = 975, a in = 40 mb for RHIC 18] and (AT coU ) = 1670, a in = 72 mb for the LHC [H^ |. 
^ilPp~ ^ s taken from Eq. (|3ufl but with the photon fragmentation function accounting for the jet energy loss. 

We will assume, as we did for pion production, that photon production via fragmentation of a jet occurs after the 
jet parton leaves the QGP. Therefore, the photon fragmentation function including the full spatial distribution and 
the secondary jets, is given by Eq. H29JI . with the substitution ir° — ► 7. 

The preequilibrium contribution of photons, corresponding to photons emitted after the transit time of the two 
nuclei but before thermalization time, is not explicitly included in this work. However, a rough estimate might be 
had by choosin g a smaller formation time, as we do. The modeling of those contributions is accessible by the parton 
cascade model 55] . Finally, in order to have a complete photon description, we have also calculated the background 
production, which mainly comes from the decay ir° — > 77. This is given by [56j 

dN^BG _ f ^0 . dN*° dP{p v 



- dy dp*. Q . (53) 

dyd 2 pr J dy lr °d 2 plf dyd 2 p T 



All the previous procedures for jet energy loss, initial spatial distribution and the effect of secondary jets are included 
in t 
by 



in the calculation of the pion spectrum — =£— — n - - In the pion center of mass frame, the photon distribution is given 

dyT d 2 pZ, 
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where 

p7 , / . 2 a , (£> cos - |p^o|) 2 

£J m = Pt cosh / sin 6» H 5 ( 55 ) 

V m 7T° 

and is the angle between p^o and p 7 . With the rj branching ratio r^ 77 /T n ~ 40% [FT) and its relative yield 
TV 1 ' /TV 71 " ~ 0.5 [58|. we simply multiply Eq. (|53() by a factor of 1.2 to include the rj contribution. 



E. Results 



Each contribution of high px photons, except the background, are shown in Fig. ^| for central collisions at RHIC 
and the LHC. The energy loss is included in all processes involving jets. Prompt photons have been split into the direct 
component (N-N) and the fragmentation component. All those processes, except the contribution from jet-medium 
bremsstrahlung, have been presented in Ref. |39j for the case of no energy loss. For RHIC, as in Rcf. 39], the high 
Pt region is dominated by direct photons. However, in Ref. |39j . the jet-photon conversion was dominant below 6 
GeV, while in our study, direct photons dominate all the high pt spectrum. A few factors explain this difference. The 
jet e nerg y loss is included here; the constants C ann and Cc om appearing in Eq. (|37f) have been set equal to 1.916 in 
Ref. [33 ; the Kj et factor in the original publication is larger than ours: Kj et = 2.5 is used for both RHIC and the LHC 
while we use Kj et = 1.7 for RHIC and 1.6 for LHC. Finally, no K 1 factor has been used for the direct contribution in 
Ref. [3{j. It is however satisfying that the inclusion of jet energy loss does not spoil the original premise: jet-photon 
conversion is an important source of electromagnetic radiation. 

At the LHC, our result is dominated by direct photons for pt above 20 GeV, but there is a window, below 14 
GeV, where the jet-photon conversion in the plasma is the dominant mechanism of photon production. In Ref. |39j |. 
however, the jet fragmentation (called bremsstrahlung in their study) was the most important process at the LHC, 
but jet suppression was not included. Photon production via jet bremsstrahlung in the plasma (dotted lines) turns out 
to be weak, but non-negligible. It is approximately a factor 3 below the jet-photon conversion contribution. Finally, 
the thermal induced photons (short dashed lines) are far below all other contributions in intensity. 

New results for total photon production, after background subtraction, are now available |59| . Our calculations are 
compared to experimental data on the left-hand side of Fig. ^] The solid line includes the prompt photon contribution 
(pQCD), the QGP (jet-th, th-th and jet-bremss.) and hadron gas contribution, extracted from Ref. 0. This latter 
contribution becomes important only for px < 1.5 GeV. The initial condition for the thermal phase corresponds to 
Ti — 370 MeV and Ti=0.26 fm/c. For a better comparison with data, we have extended our calculation down to 
Pt=1 GeV. No cutoff has been applied on either jet-th or pQCD process. NLO calculations are not very reliable is 
this region, but thermal induced reactions and hadron gas contributions turn out to dominate here and NLO results 
play a minor role. When the jet-photon conversion is not included (dot-dashed line), the total photon production is 
reduced by up to 45 %, around pr — 3 GeV, showing the importance of that process. The result expected from N-N 
collisions scaled to Au-Au is also shown (dashed line). The plasma contribution, specially the jet-thermal process, is 
very important for pt < 6 GeV. However, large error bars prevent a strong claim about the presence of a QGP. 

Fig. 1151 shows the QGP photons for three different initial conditions: (T; = 447 MeV, 71=0.147 fm/c), (7* = 370 
MeV, Ti=0.26 fm/c) and (Ti — 280 MeV, Ti=0.6 fm/c). As the high-p^ photons are produced early in the collision, 
they may be affected by the initial conditions. However, the high-py region is dominated by jet-therm processes, 
which are weakly sensitive to (Ti, Ti), since the jet distribution function fj et has a weak temperature dependence. We 
see that going from Ti=0.6 fm/c to Ti=0.26 fm/c increases the photon production in the QGP phase by less that a 
factor 2; this additional contribution could be interpreted in some sense as a preequilibrium contribution. 

Our prediction for the LHC is shown in the right-hand side of Fig. The signature of the QGP phase is much 
stronger than at RHIC, increasing the photon yield, relatively to N-N scaled to A- A collision, by one order of magnitude 
around pt = 3 GeV, where the hadron gas contribution turns out to be negligible. 

Finally, we have calculated the ratio of the total number of photons and the background photons 



-j2 — 2 r S a ^ other sources 

7Total/7BG = ^ ITJVT^ ( 56 ) 

d' 2 p T dy 



and compared in Fig.^J with the result from PHENIX [6(j, with and without the QGP contribution. The calculation 
including the QGP contribution is in agreement with the data from RHIC, except for a few data points in the range 
7 < pt < 9 GeV. Without the thermal contributions, the resulting line (dot-dashed) does not overlap at all with 
the experimental data. That could constitute a signature of the importance of the jet-photon conversion inside the 
QGP, since this is the most important thermal process as we have seen in Fig. ^| We also show the weak effect of 
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FIG. 13: (Color online) Contributing sources of high-pT photons at mid-rapidity in central Au-Au collisions at RHIC (left panel) 
and Pb-Pb collisions at the LHC (right panel). Solid line: jet-photon conversion in the plasma; dotted line: bremsstrahlung 
from jets in the plasma; short dashed line: thermal induced production of photons; long dashed line: fragmentation of jets 
outside the plasma; and dot-dashed line: direct contribution from the primordial hard scattering. 
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FIG. 14: (Color online) Total production of photons in central Au-Au collisions at RHIC and Pb-Pb at the LHC. The solid lines 
include all process from Fig. 1131 and the hadron gas contribution 1411 . while the dot-dashed lines do not include the jet-thermal 
contribution, p-p collisions scaled to A-A are shown by the dashed lines. Data at RHIC are from PHENIX |5fl^ . 



the initial temperature. The ratio 7Totai/7BG at T = 447 MeV is only enhanced by ~ 5% relatively to the result at 
Ti=370 MeV. Finally at the LHC (right panel), the thermal contribution is also visible: including the photons from 
the thermal phase enhances the calculation by ~ 15%. 



VI. SUMMARY AND CONCLUSIONS 



We have used a complete leading-order treatment of jet energy loss in the QCD plasma to calculate the pion and 
photon spectra for both RHIC and the LHC. The calculations have been confronted with available data from RHIC 
and turn out to be in good agreement. These results reinforce the idea that high px suppression is a final state effect 
caused by jet energy loss through bremsstrahlung in the hot medium. 

The neutral pion nuclear modification factor at RHIC as been reproduced with an initial temperature Ti =370 
MeV and a formation time 7^=0.26 fm/c, corresponding to dN/dy=1260. Those parameters are consistent with the 
analysis in |4l|. Raa has shown a large dN/dy dependence, but a weak dependence on the initial temperature Ti 
(provided the starting time Ti is changed to keep dN/dy constant). The calculation included the nuclear geometry; 
assuming that all jets are produced at the center overestimates the suppression by ~ 50%. 

We have also computed the production of high pr photons from the initial collision, from the medium, and from 
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FIG. 15: (Color online) Production of photons during the QGP phase for three QGP initial conditions: solid line, (T; = 447 
MeV,Ti=0.147 fm/c); dashed line, (Ti = 370 MeV,r;=0.26 fm/c); and dot-dashed line, (T; = 280 MeV,-n=0.6 fm/c). 
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FIG. 16: (Color online) Ratio of all photons over all decay photons as a function of pr, in Au-Au collisions at RHIC (left panel) 
and Pb-Pb at the LHC (right panel), with and without the thermal contribution. In the left panel, the effect of the initial 
temperature is also shown. The data are from the PHENIX [6(J collaboration. 



jet-medium interactions. The jet-medium photons improve the agreement between experiment and theory at RHIC, 
and they are expected to dominate the signal at the LHC below about 14 GeV. Thermal photons from the medium 
are not very important to either experiment, in the kinematical range on which we have concentrated. In light of 
these results, the in-medium production of dileptons should also be reconsidered. Work on this topic, and others, is 
in progress. 
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